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We study the thermodynamic condensation of microcavity polaritons using a realistic model of 
disorder in semiconductor quantum wells. This approach correctly describes the polariton inhomo- 
geneous broadening in the low density limit, and treats scattering by disorder to all orders in the 
condensed regime. While the weak disorder changes the thermodynamic properties of the transition 
little, the effects of disorder in the condensed state are prominent in the excitations and can be seen 
in resonant Rayleigh scattering. 
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Considerable effort has been recently devoted to the 
realisation of a Bose-Einstein condensate of polaritons 
in III-V and II- VI semiconductor microcavities [3, 0, . 
The very light mass of these composite bosonic particles 
promises relatively high transition temperatures, estab- 
lishing these systems as ideal candidates for observing 
condensation. A significant challenge to the realisation 
of an equilibrium condensate might be represented by 
the short polariton lifetime (caused by the finite quality 
of cavity mirrors) and by the suppression of thermaliza- 
tion processes by acoustic phonons at small momenta - 
the 'bottleneck effect'. However, very recent develop- 
ments have suggested that, by positively detuning the 
cavity energy above the exciton energy, and by increas- 
ing the non-resonant pump power to increase particle- 
particle scattering, thermalization can be dramatically 
amplified While unambiguous evidence for equilib- 
rium condensation still remains uncertain, much progress 
has been achieved in this direction. This includes the 
observation of a non-linear threshold behaviour in the 
emission intensity under non-resonant pump, the de- 
crease, above threshold, of the second order coherence 
function 2] together with a characteristic change in the 
momentum space distribution and, recently, interference 
patterns in far-field emission have been measured 0|. 

Theoretical effort has also been devoted to predicting 
properties and possible signatures of polariton conden- 
sation In this Letter we consider how disorder, 
through the distribution of excitonic energies and oscil- 
lator strengths, affects such signatures. Even with the 
sophisticated growth technologies used in current struc- 
tures, the presence of interface and alloy disorder in- 
duces noticeable effects. Here, we will show that the 
response under resonant Rayleigh scattering (RRS), the 
coherent scattering by the disorder of an injected pho- 
ton into directions other than its original direction, pro- 
vides a unique probe of the condensed regime. Already 
in the low density (linear) regime, disorder determines 
the RRS response 0, and the inhomogeneous broad- 
ening of the polariton photoluminescence (PL). Using 



a quantitatively accurate model for the exciton disor- 
der |3, |l£| that has already been well-tested in the linear 
regime of excitation, we investigate the effects of disor- 
der on an equilibrium polariton condensation. We find 
that the thermodynamic properties of the transition are 
weakly affected by small disorder, as expected. However, 
at densities above the threshold expected for condensa- 
tion, the normal modes supported by the cavity change 
from the lower and upper polariton modes to new col- 
lective excitations Accordingly, the response of the 
condensate to an additional RRS probe also changes. In 
noticeable contrast with the non-condensed regime, the 
frequency-resolved RRS emission from the linear Gold- 
stone mode exists both above and below the chemical 
potential. Moreover, the spectrum exhibits features di- 
rectly related to the disordered quasi-particle spectrum. 
Here, in analogy with the BCS theory, the quasi-particles 
are given by 'particle-hole' excitations (i.e. bound exci- 
tons) coupled to the coherent photon field via the ran- 
dom, disorder dependent, oscillator strength. 

The linear response of a resonantly pumped polaritonic 
system to an external disorder potential has been recently 
studied in [ll[. In that work, in contrast with the case 
analysed here, the coherence of the system is driven by 
the pump, and moreover the effect of disorder is included 
at a perturbative level. 

The influence of quantum well disorder on excitonic 
energies and oscillator strengths has been much studied 
in the last two decades (for an exhaustive discussion see, 
e.g., 9] and Refs. therein). Here, we assume the dis- 
order to be correlated on a length scale £ c larger than 
the exciton Bohr radius a x = e/e 2 /i r , where \i r is the re- 
duced mass (henceforth we will set h = 1). Accordingly, 
we factorise the in-plane relative and centre of mass co- 
ordinates, \P a (r e , rft) ~ y?i s (r)$ a (R), and focus on the 
excitonic centre of mass motion, 
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where the energies are measured w.r.t. the exciton band 



edge i.e., the band gap minus the exciton Rydberg 
1Zy x = (2fi r a 2 )~ 1 . The effective disorder potential V(R) 
represents the microscopic structural disorder averaged 
over the electron-hole motion 0] and can be approx- 
imated, e.g., with a Gaussian noise, (V(R)V(R/)) = 
(cr 2 £ 2 c /L 2 ) Y}^ c e ic ^ R - R/ \ where L 2 is the quantisation 
area. 

Even though, in two dimensional non-interacting sys- 
tems, all states are localised by the disorder potential, 
the character of the excitonic wavefunction changes sig- 
nificantly from below to above the band edge E x . The 
coupling strength of an exciton to light changes accord- 
ingly. For dipole- allowed transitions, the exciton oscilla- 
tor strength g a ^ p is proportional to the probability am- 
plitude (pis(0)®a,p of finding an electron and a hole at 
the same position and with centre of mass momentum 
equal to the photon momentum p 0, , 
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where is the dipole matrix element. Here, we solve 
Eq. numerically on a grid of 120 x 120 points (for 
which convergence is reached) for a system of size L = 
1/im, a = 2meV and £ c — 166 A. From the evaluated 
eigenvalues e a and eigenstates $ a ,p, we can derive the 
excitonic density of states DoS(e), the coupling strength 
g ajP and its squared average (Fig. [I}: 
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where (...) is the average over different disorder reali- 
sations. This quantity is related to the excitonic optical 
density by D(e) = DoS(e)g 2 (£,0). 

We now consider the following Hamiltonian describing 
excitons with random energies e a dipole coupled via g ajP 
to the cavity field i(j p : 



~jjlJ2J2 (9a iP ^pbia a +h.c.) . 



(4) 



Interactions are approximated by exclusion, so each ex- 
citon level e a is modelled by an electron-hole pair a 
and 6, with the total occupation restricted to one, i.e. 
b^b a + a^ a a a = 1. For the thermodynamical properties 
of this model (e.g., the critical temperature), this is a 
good assumption at low enough densities, where only the 
strongly localised (Lifshitz) states in the tail are pop- 
ulated. The dispersion for photons in a microcavity 
of width L w is approximated as parabolic, uo p ~ uoq + 
p 2 /2m p h, where c^o = lucj L w yfe and m p h = 2i\yJ~ej cL w . 
It is convenient to rescale the exciton-light coupling g a ^ p 
according to g ajP i— > g aj p\/^Y^ ^x/27r, where N = 




S 

-2 



e=p /2m 
g 2 (£,p) 





1 








pi- 









6 

P/Pstep 



10 



12 



FIG. 1: (a) Plot of the squared coupling strength |g«,p| 2 ver- 
sus energy (160 realisations of disorder) for p = 0; (b) con- 
tourplot of the squared averaged oscillator strength g 2 (e, |p|) 
versus energy and momentum. The resolution in momentum, 
Pstep = 2ty/L = 6.3 x 10 4 cm -1 corresponds, for a cavity of 
ujo — 1.68eV, to an angle of # s te P = tan -1 (cp s te P /^o) = 36°. 
The free particle dispersion |p| 2 /2ra x (solid line), and the 
squared averaged oscillator strength for two values of mo- 
menta, |p| = and |p| = lOpstep [0 — 82°) (plus symbols) are 
explicitly plotted. 



7£y x L 2 ra x /27r is the inverse level spacing measured in 
units of the excitonic Rydberg energy. This corresponds 
to measuring the density of particles in units of the Bohr 
radius squared. In thermal equilibrium, the total num- 
ber of excitations, TV = J2 a ( b i b a + a « a l)/ 2 + J2 P ^I^p, 
can be fixed by introducing a chemical potential, \i. 
The dimensionless density p = (N)/N corresponds to 
((7V)/L 2 )a 2 47r/i r /ra x , where {N)/L 2 is the density of 
particles per area. 

Making use of standard path integral techniques, the 
mean-field equations for the static and uniform photon 
field x/j and for the total number of excitations can be 
simultaneously solved in a similar way as in Refs. 
However, in contrast to previous work, we introduce here 
a realistic description for the excitons in the disordered 
quantum wells and make use of the energies and coupling 
strengths evaluated numerically. In addition, the aver- 
age over different disorder realisations outside the energy 
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FIG. 2: Mean field phase diagram for the dimensionless crit- 
ical temperature 2T c £;b/Qr versus the dimensionless den- 
sity p for effective zero detuning c^o — E x = — 0.94meV and 
Q K = 26meV. A detail of the low density region is shown in 
the inset (a) , while the plot of the mean- field order parameter 
i/j/y/N versus the density for Tk B = 20K (2Tk B /^R = 0.13) 
is shown in the inset (b). 



interval evaluated numerically (Fig.^a)) is taken by ex- 
trapolating the numerics: in the low energy Lifshitz tail, 
we approximate the distribution of \g a ,o\ 2 with a delta 
function at its extrapolated mean value g 2 (£,0), while, 
in the high energy region, we use the Porter-Thomas dis- 
tribution V(x = |^q;,o| 2 ) = exp[— xj (2x)\j \j2nxx, where 
x = g 2 (e,0). The scale of |^ a ,o| 2 is set so as to fix the 
integrated optical density J deD(e) = Or/4, where Or 
is the measured Rabi splitting. In the resulting mean- 
field phase diagram (Fig. [2j, increasing the density of 
particles at a given temperature, the system goes under 
a phase transition from a non-condensed phase (ip = 0) 
to a phase where the polaritons condense in the lowest 
momentum state (-0^0). 

The incoherent PL spectrum and the RRS response 
are evaluated by considering fluctuations of the photonic 
field above the mean-field solution: 
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Here, when Uh = ^irksTh ^ 0, the matrix elements of 
the inverse quasi-particle Green's function are given by 

K pI (vh) = £ P ,q fab. + a) P ) + ^ 9*, P 9*,<i (6) 
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where h{x) = t&nh(x /k B T)/[x(ujl + 4x 2 )], E a = 
y/(s a /2) 2 + \g ai o\ 2 i/j 2 /N is the energy of an exciton in 
a coherent field, e a — e a — n and uj p = uo p — fi. The 
quasi-particle Green's function can be decomposed into 
(momentum) diagonal and off-diagonal contributions, 
(1 - 2 ^a _l K p l ^° . The off-diagonal terms, 
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as they break translational invariance and therefore have 
a zero average over different disorder realisations, can 
be neglected when evaluating the incoherent PL inten- 
sity P(cj,p) = tib((^)W(uj,p), where tib(w) is the 
Bose occupation factor and 
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is the spectral weight. However, allowing the normal- 
modes supported by the cavity to scatter via their exci- 
tonic component in directions different from the incoming 
one, these terms are essential in finding the RRS inten- 
sity: 
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Here, the filter function F p = \\K^\ 2 - [K^] 2 \~ 2 de- 
scribes the propagation of the injected and detected pho- 
tons via the normal modes supported by the cavity, while 
the scattering function S pci 
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describes the probability to scatter, via the excitonic 
component, from the injected momentum p to the de- 
tected one q. The incoherent PL P(u,p) and spec- 
tral weight W(w,p), as also the disorder averaged RRS 
response (J pq (a;)) are shown in Fig. [3| for two values 
of the density at a fixed temperature, showing both 
the non-condensed (Figs. Efa,c,e)) and condensed case 
(Figs.Etb,d,f)). 

The RRS signal reflects many of the features already 
present in the spectral weight. The similarity is due to 
the filter functions, which limit the RRS response to the 
normal modes supported by the cavity and are respon- 
sible for the ring-shaped emission observed in experi- 
ments 0, 0. When uncondensed 0), = 0, 
the filter terms \K^ d \~ 2 coincide with the squared po- 
lariton Green's function and the scattering term lii^ 1 ) ! 2 
gives the probability for an exciton to scatter from p to 
q. At ultralow densities, the model used here is equiv- 
alent to that used in Ref. and the RRS response 
is in agreement with that of Refs. 0, However, 
when condensed, the polariton modes are replaced by 
new collective modes, the lower polariton becomes a lin- 
ear Goldstone mode, and two new branches appear be- 
low the chemical potential, which are seen as gain in the 



4 



Non-Condensed 



Condensed 




0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 



FIG. 3: Contourplot of the incoherent PL P(u;,p) (a,b), the 
spectral weight W(w, p) (c,d), and the disorder averaged RRS 
intensity (ipq (<<;)) for |p| = |q| (e,f) versus the rescaled en- 
ergy 2 (u) — /jl)/Qr and the dimensionless momentum |p|a x , 
for cj - E x = -0.94meV, Q R 26meV, fc B T 20K and 
= 0.01: non-condensed (left column) (p ~ 0, /i = 
— 16.5meV) and condensed (right column) (p = 3.6 x 10 -3 , 
p=-11.3meV). 



spectral weight. These changes are also seen in the RRS 
response, including RRS response at energies below the 
chemical potential (Figs. E{f)). 

These noticeable changes could be observed in both PL 
emission, the product of the spectral weight times the 
Bose occupation factor ^(w), and the RRS response. 
However, because of the occupation factor, PL emission 
from above the chemical potential is exponentially sup- 
pressed. In addition, the PL as it is plotted in Fig. [3| 
excludes the emission from the condensate, which in ex- 
periments might obscure these features. In contrast, RRS 
response is not weighted by occupation, but instead by 
the scattering function, .S'pq, which depends on the mod- 
ulus squared of the Green's function for an exciton in 
a coherent field. As a consequence, RRS represents a 
unique probe for observing condensation in polariton mi- 
cro cavities. 

Using the full distribution of oscillator strengths has 
observable effects in the condensed state and is vital in 
our treatment. When uncondensed, while the spectral 
weight depends only on the excitonic optical density (and 
thus on the average of |^ a ,o| 2 ) 5 the RRS depends also on 
Ig^ol 4 . This determines a sharper energy dependence 



in RRS than in optical density. However, when con- 
densed, both optical responses are determined by the the 
full distribution of oscillator strengths, rather than only 
the mean squared oscillator strength and its mean fourth 
power. This is because the energy E a of an exciton in 
the presence of a coherent field leads to a dependence 
on the full distribution of |#a,o|- The density of states 
of these excitonic quasiparticles can be directly seen in 
the spectral weight between the upper and lower polari- 
ton modes. With a constant g, the minimum of E a is 
sharp, occurring for excitons where e a = 0. In contrast, 
with the distribution of \g a ,o\ used here, the minimum of 
E a occurs for states where there is first a finite proba- 
bility of having negligible coupling to light. Below this 
energy, other than the linear mode, there is a suppres- 
sion of the spectral weight. The energy dependence of 
the distribution of coupling strengths is also important, 
and is responsible for the decrease of polariton splitting 
at higher densities. As the chemical potential increases, 
the excitons with energies close to the chemical potential 
couple less strongly to light, and the effective splitting 
decreases. 

The treatment in this Letter considers only disorder 
acting on the excitons; it has been suggested that 
disorder of the mirrors, acting on the photons might also 
be relevant. The inclusion of such disorder can be ex- 
pected to modify the scattering term, but the filter func- 
tions are expected to remain unaffected. 

To conclude, we have described how polariton conden- 
sation in semiconductor microcavities can be investigated 
by RRS. By making use of a realistic model of disorder 
in quantum wells, we identified the changes which can 
be observed in the RRS spectrum when, with increasing 
density, the system, still in the strong coupling regime, 
crosses the the phase boundary to a condensed state. Be- 
cause RRS, in contrast to PL, is not weighted by the Bose 
factor, and because RRS does not contain a strong sig- 
nal at the condensate frequency, which could obscure the 
subtle features in PL, RRS provides a promising probe of 
polariton condensation in semiconductor microcavities. 
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